T-violation in flavour oscillations as a test for relativity principles at a neutrino factory 
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We study the effects of violation of tfie equivalence principle (VEP) or violation of Lorentz invari- 
ance (LIV) in the neutrino sector on the asymmetry between T-conjugate oscillation probabilities, 
APt = Piva vp) — Piyp Va), iu a three-flavour framework. We find that additional mixing 
due to these mechanisms, while obeying all present bounds, can lead to an observable enhancement, 
suppression, and/or sign change in APt for the preferred energies and baselines of a neutrino fac- 
tory. The measurement of this asymmetry can be used to establish a new upper limit of order 10~^® 
on VEP or LIV in the (//e, i^^) and (j/o, j/^) sectors. 
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I. REVIEW AND MOTIVATIONS 

The phenomenon of neutrino flavour oscillations fol- 
lows from the mixing of nondegenerate neutrino states. 
In the "standard" scenario, a flavour neutrino is a 
linear combination of neutrino mass eigenstates and 
the degeneracy is broken by means of small neutrino 
masses. Alternative mechanisms that do not require 
neutrino masses (or more precisely, nondegenerate neu- 
trino masses) abound. One possibility is through a vio- 
lation of the equivalence principle (VEP) in the neutrino 
sector 0, 0]. Here, a flavour neutrino is a mixture of 
neutrino gravitational eigenstates and flavour-dependent 
neutrino-gravity couplings serve to lift the degeneracy. 

An essential phenomenological difference between the 
standard and the VEP mechanisms is the manner in 
which the oscillation probability P{i'a ~^ vp) depends 
on the neutrino energy E. In the two-flavour CcLS6, cl for- 
mal transformation from the former to the latter scenario 
is readily accomplished by replacing in P{va vp): 



to the LSND 



-18 



5w? 
'2E 



(1) 



where Svn? = rn^—m\ is the neutrino squared mass differ- 
ence, (57 = 72 — 7i the difference in the neutrino-gravity 
couplings, and (p the Newtonian gravitational potential. 
Note that a violation of Lorentz invariance (LIV) , where 
the breaking of the neutrino degeneracy is assumed by 
dissimilar neutrino asymptotic speeds Vi 's^, may also 
contribute to the oscillation phenomenon with the same 
energy dependence as in the VEP case, save for a minor 
notational change: \(l)\6^ — > 5v/2, where 5v = V2 — wi ^^j. 
Henceforth, wc shall use only the VEP convention. 

To date, the most comprehensive set of bounds on the 
VEP parameter \(t)\5^ comes from the CCFR experiment 
0: \(t>5^\ ^ 10~^^ for all three oscillation channels in 
the two-flavour limit 0. These extremely severe con- 
straints have since excluded VEP-LIV as an explanation 



transition puzzle which re- 
quires \(t)5^\ ^ 5x 10~^® |3, as well as Mikheycv-Smirnov- 
Wolfenstein (MSW) type solutions (with \(l)5-i\ > 10"^^ 
0,0,01) to the solar neutrino problem (see later). 

In the case of atmospheric i/f^ and 17^ disappearance, 
the totality of the Super-Kamiokande atmospheric neu- 
trino data im strongly favours an energy dependence 
of 1 / E" , where n > 0, in the oscillation frequency 
Hi nil. No acceptable VEP-LIV solution exists for 
the complete set of contained, partially contained, and 
upward going muon events. In an interesting analysis, the 
authors of Ref. considered VEP-LIV as a sub-leading 
effect to standard mass- induced oscillations within a two- 
flavour ^ i>T framework, and derived the upper limit 
\(j)Sj\ < 3 X 10-24 at 90 % CL. for unconstrained VEP 
mixing in the {v^, Vt) sector. For maximal VEP mixing, 
the bomid on \(t>5^\ can be as stringent as ^ lO^^^. 

On the other hand, it has been demonstrated in a num- 
ber of recent works jlj, [la, that vacuum oscillations 
induced by a "just-so" VEP with \4)6^\ ~ 10"^'* and 
sin^ 29 ~ \ can still provide a good description of the 
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solar neutrino data ^ M, M, 2^ 2^\M\M- The anal- 
yses of Ref. put the solution's goodness of fit at 60% 
C.L., a figure that is comparable even to the confidence 
level of 75 % for the standard large mixing angle (LMA) 
solution. However, the LMA oscillation parameters have 
now been confirmed by the terrestrial experiment Kam- 
LAND, through the observation of reactor Vf. disappear- 
ance over a typical baseline of 180 km While these 
new results do not exclude VEP at the aforementioned 
scale, a pure VEP-LIV solution to the solar deficit is 
no longer justified. 

It is clear from the above discussion that neutrino os- 
cillations provide very stringent tests for the principle 
of equivalence and Lorentz invariance. We examine in 
the present work the possibility of further constraining 
VEP-LIV in the neutrino sector with future experiments. 
Currently operating long baseline terrestrial experiments 
(e.g., K2K), those under construction (e.g., MINOS), and 
prospective neutrino factories psl ] will provide for us a 
means to obtain more precise information on the neutrino 
mixing properties, as well as an opportunity to establish 
CP- or T-violation in the lepton sector. In particular, 
the capacity of a neutrino factory to produce both elec- 
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tron and muon neutrino and antineutrino beams in high 
intensities means that it is in principle possible to di- 
rectly compare the T-conjugate oscillation probabilities 
Pii^fi and P{ve Vfj) '2^. Despite the associated 

experimental challenges, this is theoretically the clean- 
est way to resolving any CP- or T- violating phases in 
the neutrino mixing matrix: a nonzero measurement of 
APy = P{i'a i^fs) — P{^i3 ~^ i^a) ncccssarily implies 
the existence of such a phase, in contrast with the CP- 
asymmetry APpp = P{i'a ~^ i^f}) ~ P{pa ~^ ^js) which 
inevitably contains additional contribution from terres- 
trial matter effects.^ Finally, a neutrino factory is ideal 
for this study since it involves high energy neutrinos trav- 
elling a very lon g di stance, which tends to enhance the 
VEP-LIV effectsji^. 

Assuming three neutrino flavours and that standard 
oscillations account for both the atmospheric neutrino 
anomaly and solar neutrino deficit, we consider VEP- 
LIV as a sub-leading mechanism (while respecting all 
present bounds), and examine its consequence on the T- 
asymmetry APj^ for the design parameters of a neutrino 
factory. We find that although the presently allowed 
|(/)|(57's are far too weak to produce any effect on the ac- 
tual oscillation probabilities beyond the few percent level, 
the T-asymmetry can be significantly enhanced, sup- 
pressed, or flip sign, depending on the VEP-LIV breaking 
scale as well as the energy of the neutrinos. 

We shall begin by introducing the parameters relevant 
for the present phenomenological study 

II. DEFINITIONS OF PHENOMENOLOGICAL 
PARAMETERS 

Consider the evolution of a three-flavour system in a 
medium, 

= (2) 

where v = (z^e, v^, VrY^ . Suppose that, in addition 
to the usual mass- mixing and matter potential terms, 
UmMU\.^ and V , the total Hamiltonian H governing this 
evolution contains also a third term UgGUq arising from 
a violation of the equivalence principle (VEP): 

H = ^UmMUI + 2E UgGuI + V, (3) 

where 

M ~ Diag(m^, TOj, Tig), 
G = Diag(7i, 72, 73), 

V = V2 GpNe X Diag(l, 0, 0). (4) 



Here, 2 3 a-re the squared masses of the neutrino mass 
eigenstates, (j) is the gravitational potential, 71,2.3 are the 
neutrino-gravity couplings, E is the neutrino energy, Gp 
is the Fermi constant, is the electron number density 
of the medium, and Um,g are two independent unitary 
matrices relating the flavour eigenstates to, respectively, 
the mass and gravitational eigenstates. Note that the 
subscript capital M stands for "mass" , and is not to be 
confused with the commonly used small m which denotes 
"matter" . 

A generic 3x3 unimodular unitary matrix may be 
parameterised with three mixing angles and six phases: 

U = e"' Diag(l, e*", e"^) U' Diag(l, e'«, e'^), (5) 

where 77, a, ^, and x are five arbitrary phases, and 
the matrix [/' contains the three mixing angles and the 
sixth phase, i.e., the CP- or T- violating phase, such as 
the familiar CKM matrix.^ The phase is an unobserv- 
able overall phase. ^ and x may correspond to phys- 
ical Majorana phases (if neutrinos are Majorana par- 
ticles), or to unphysical phases (if neutrinos are Dirac 
particles), neither of which contribute to flavour oscilla- 
tions. If only the mass term is present, the remaining 
two phases a and /3 are unphysical, and can always be 
absorbed into the deflnitions of the flavour eigenstates, 
e.g., I'^'j = e"'^"^' ffj^, leaving U'j^j the only observable en- 
tity. Similarly, if UgGUq is the unique term, only Uq 
contributes to the oscillation phenomenon. In the gen- 
eral mass plus gravity case, however, it is not possible to 
completely eliminate the a and /3 phases. If we choose 
Um = U'j^j, it follows that the total Hamiltonian H must 
take the form 

H = -^U'mMU^i + WUaGU'}w'^ + V 
2E 

= H^' + WH^W^ + V 

= H + WH^W\ (6) 

where W = Diag(l, e'", e*'^), with a = ac — cum and 
fi = j3G — Pm, and we have defined H = H^^ + V as 
the "standard" mass-induced mixing plus matter poten- 
tial sub-Hamiltonian for future reference. To transform 
from to WH^W\ one simply makes the following 
replacements: 







- 














H-G p«(/3-") 




ttG t(a^l3) 



Thus the inclusion of a gravitational term in the total 
Hamiltonian leads to, in principle, eight, not six, extra 
observable parameters! 



^ It should be noted that T- violating effects can arise from an earth 
matter distribution that is not symmetric about the midpoint of 
the neutrino's traj ect ory, indepe nde ntly of the intrinsic CP- or 
T- violating phase l27l . See Ref. |28|| for a recent study. 



^ A common practice is to express some of these extra phases in 
terms of SU (3) generators, utilising the diagonal Gell-Mann ma- 
trices A3 and Ag. We find our simple parameterisation to be 
more convenient for the present phenomenological study. 
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In the following, we shall make the simplifying assump- 
tion that the splitting between two of the three gravita- 
tional eigenstates is too small to be probed by present day 
experiments and experiments in the near future. This 
near degeneracy means that only one \(l)\5^ and two VEP 
mixing angles contribute to the observable part of H'^ , 
and the CP- or T- violating phase inherent in [/^ can be 
absorbed into the definitions of the phases a and f}. We 
shall adopt the convention 71 ~ 72 and ^7 = 73 — 72 , and 
parameterise the effective U'q as 

cos If sin If \ 

U'n — I — sin sin cos?/' sin ^ cos (/3 , (8) 
— cos ■(/' sin (f — sin -0 cos '0 cos f I 

where f> and "0 run from — 7r/2 to 7r/2. The resulting 
three-flavour sub-Hamiltonian H"^ takes the form 

2sin^(/3 sin'0sin2iy9 cosV'sin2(/? 
-E'|(/)|(57 I sin'0sin2iy9 2 cos^(/?sin^^ cos^</3sin2'!/) 
cosV'sin2(p cos^((9sin2-0 2 cos^t/J cos^tA 

(9) 

which reduces to a pure two-flavour description in terms 
of one mixing angle when either or is fixed at 0, 7r/2, 
or — 7r/2. Specifically 

svcP'i}) ~ \ i^e ^ VEP mixing, 
sin^^A = ^ Ve Vt VEP mixing, 
sva^f = ^ ^ Vt VEP mixing. (10) 

Observe that the sign of the parameter \(l)\5^ may be 
positive or negative. As we shall see, this has important 
consequences when considered together with H^^ -t- V . 

For the mass- mixing term , we adopt the usual 
MNS parameterisation [sOl for U'j^: 



C1C3 
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S1C3 
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C1C2 - siS2S3e S2C3 

-C1S2 - SiC2S3e*'^ C2C3 



S1S2 - CiC2S3e 

(11) 

where Si = shiOi, Ci = cos9i, i = 1,2,3, and 6 is the 
CP- or T-violating phase. The solar and atmospheric 
oscillation parameters are identified as 



Sm, 



atm 
2 

sun 



'32 



mo 



r 2 _ 2 
Om2x = Tlj 



'3 ~ ^2 J ''atm = [72 , 

ll, ^'sun = 01, (12) 

assuming that maximal <-^- oscillations with 
(Jm^tm — 3x 10"'^ eV^, and the large mixing angle (LMA) 

,2 ^, c ^, 1 n — 5 n-^fS 



solution {Sm^^^ ~ 5 x 10"^ eV^, sin20sun — 0.8) explain, 
respectively, the atmospheric and solar neutrino anoma- 
lies. Given these solutions, a three-flavour analysis of the 
CHOOZ data sees the remai ning angle 6*3 constrained to 
tan2 03 < 0.055 at 90% C.L. 



where A 1,2, 3 are the eigenvalues, and the energy depen- 
dent unitary matrix U contains by definition one in- 
dependent effective CP- or T-violating phase. Assum- 
ing constant Ng and 0, the asymmetry between the 
T-conjugate oscillation probabilities P(;/q — > Vfs) and 
P{vf3 Va) is given by 

APt = P{va vp) - P{v0 va) 

J . Ai2i . A23X . A31L 

— IbJ sm sm sm , 14) 

2 2 2 ^ ' 

where = — A^, with i,j = 1,2,3, is the difference 
between the zth and the jth eigenvalues, and 



J = Jjp = Im{Ui3iU^jU*iUaj) 



(15) 



is the effective Jarlskog factor |33 encapsulating the said 
effective phase. 

Explicit evaluation of the quantities J and Ay is a 
tedious, though straightforward, task. Fortunately, there 
exists a very useful identity which relates them to the 
original Hamiltonian in flavour space js^ : 



A12A23A31 J = lm{He^H^rHre)- 



(16) 



Observe that (i) the right hand side of Eq. (|16|l is depen- 
dent only on the off-diagonal elements of H, and (ii) the 
left hand side of Eq. Hlt)|) bears a striking resemblance to 
the right hand side of Eq. ifT^ . These immediately sug- 
gest that the observable APt is potentially very sensitive 
to the additional mixing terms arising from VEP. 

We do not consider the CP asymmetry APcp = 
P{va — ^ vp) — Pip a ^ ^13) in the present work. This 
measurable inevitably contains additional contribution 
from terrestrial matter effects, thus rendering its ana- 
lytical description much more complicated. 



IV. ANALYSIS 

We now analyse the effects of VEP on the observed 
T-asymmetry. Consider the ratio 



APt J sin^f^ sin ^f^ sin 
A^ ^ J sin^ sin^ sin^^' 



(17) 



where the quantities with a tilde are as defined in Eqs. 
(|14|l and and those without correspond to their 

"standard" H'^ — counterparts. Using the matter- 
invariants 



III. WHY T- VIOLATION EXPERIMENTS? 

Solution to the Hamiltonian (|3Jl may be obtained by 
first transforming to a new basis defined as 



C/ti7C/ = Diag(Ai, A2, A3 



(13) 



ln-v{H,^H,rHre) = iHH^lH^iH^l) , (18) 



and 



J 



1 Sml^Smlr^dml^ jj^.^ 
8E3 A12A23A31 



(19) 
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where 



J*^ = S1C1S2C2S3C3 sin(5 (20) 



is the Jaiiskog factor in vacuum, we rewrite Eq. H17|l as 

APt 8£;3Ai2A23A3iJ 



APt drnioSml^Sm^J^ 



12 '^'"'23'^" ''31"- 

'sin^f^ sin^f^ sin ^f^ ' 
sin^f^ sin^f^ sin 



(21) 



Rough estimates of the differences between the eigen- 
values Ay, ij — 12,23,31 may be extracted from the 
Hamiltonian (|3Jl and Eq. Q: 

-0{V2GFN,) + 0{2E\(l)\5-f), (22) 



2E 



A. 



23 



A31 



2E 



-0{2E\4>\5^), 



(23) 



2E 



^ -0{^/2GFNe)+0{2E\(l>\5-^). (24) 



Their standard counterparts A^ are given similarly, save 
for the absence of the 0{2E\(f)\5^) terms. 
Suppose now that 



2£;|^<57l«''^"'^' 



2E 



(25) 



or, equivalently, 



|<^.,|«2.5xl0-^(^) . (26) 
This immediately allows us to set A23 ~ A23, and thus 



and only A12 is substantially affected by . Equation 
(|29|l may be further simplified if the neutrino energy and 
the experimental baseline are such that 



sm ■ 



A12L A12L 



A12L Al2i 



sm ■ 



2 ' 



(30) 



or, roughly speaking, if the experimental baseline satisfies 
the requirement 



L < 



A12I Max(^, V^GpNe, 2E\(f>S-f\ 



(31) 



The outcome of this approximation is a T-asymmetry 
ratio that is related in a simple way to the off-diagonal 
elements of the original Hamiltonians in flavour space. 



APt 8i;3Ai2A23A3lJ lm{H,^H^rHre) 



(32) 



independently of the baseline of the experiment. This is 
one of the main results of the present work. 

On-resonance: A31 ~ 0. The near cancellation of 
OiV^GpNe) with other terms in Eq. and the con- 
dition together imply for Eq. (|^ that 



A12 



2E 



-0{2E\<j>\d'y) ~ A 



12, 



and thus 



. A12L . Ai2i 

sm ~ sm . 



(33) 



(34) 



Furthermore, if the experimental baseline and the neu- 
trino energy satisfy the condition 



A23i 



sm ■ 



A23i 



(27) 



For |i5rn32| ~ 3 x lO"'^ eV^, and a neutrino factory pro- 
ducing a neutrino beam of average energy 13 GeV, Eq. 
(|25|l is a good approximation for {(pS'^l ^ 4 x 10"^'', a 
region of parameter space well below the present CCFR 
bounds. Further approximations are in line, depending 
on whether A31 is on- or off-resonance. 

Off-resonance: A31 ^ 0. If there is no exact or near 
cancellation between 0{-\/2GFNe) and other terms in Eq. 
(|24|l . the condition H25() implies also that A31 ~ A31, and 
therefore 



A31L 



A31L 



Hence, 



APt 
APt 



8£;3Ai2A23A3iJ sinA^ 



'12'- 



^l,Sml, J'' Sin ^ 



(28) 



(29) 



L < 



IA31I Max(i^^sin26l3, 2E\(j)S-f\ 
then the approximations 



, (35) 



. Asii A31L 
sm — - — ~ — - — , 



. A3i£ A3i£ 
sm — - — ~ — - — (36) 



are also well-founded. Substituting Eqs. (|27ll . I|34|l and 
(|36|l into Eq. H17() . we find that the T-asymmetry ratio 
in the on-resonance case is equally well described by Eq. 
(|32|l . except for a different baseline requirement (|35|) . 

As it turns out, Eq. (|32|l is a generally valid approxima- 
tion for the preferred parameters of a neutrino factory: 
an average neutrino energy of about 10 to 30 GeV and 
a baseline roughly between 3000 and 6000 km, provided 
that Eq. 123) holds. 

It is useful to rewrite Eq. H32() as 



APt 

'aK 



1 + i?, 



(37) 
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: sin2(p = 


1 sinV=0 E=13GcV 
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FIG. 1: APt/APt versus \(j)\5"f for S = 13 GeV, maximal 
i/e ^ Vr VEP mixing (sin2(p = 1 and so? i^} = 0), and 
with the relative phase /3 set to zero. The standard oscil- 
lation parameters used here are the best fit values: |5m|2| ~ 
3 X 10"^ eV^ sin26l2 = 1, 5ml^ = 5 x 10"^eV^ and 
sin 2^1 = 0.8 . The four sets of points correspond to exact 
values of APt/APt from numerical calculations for different 
values of S3, 5, and baselines L. The solid lines come from 
the approximate equation 14311 : the top (bottom) set is for 
negative (positive) |c/>|(57's. 



with 



cyclic 

E 

ijk 



ImiHifH^iH^) Im(<i^'i^) 



(38) 

where Hf^ = {WH^W^)ij , the summation is taken over 
ijk — e/ir, fire, re 11, and we have used 



im(i/,^^ff;r;i/,^,) = im{H':^H\:,H%) = o, 



(39) 



since two of the three gravitational eigenstates are by 
assumption degenerate. The five extra parameters from 
H'^ offer a plethora of possibilities. We provide here a 
few illustrative cases, but the list is by no means exhaus- 
tive. 



A. Pure Ve i't VEP mixing 

The setting of = corresponds to pure ^ Vt 
mixing in the VEP sector. In Eq. (|38|l . the only surviving 
term is 



(40) 



FIG. 2: Same as Fig. [T] except for _E = 25 GeV. 



Using the parameterisation introduced in Sec. m we ob- 
tain 



\(j>\5j sin2(y9. 



where 



B ■ 



S2 sin(/3 — S) 



SiCi 



sin (5 



"■32 



C2 Sm^i sin/3 
S3 6111,^2 sin (5 ' 



(41) 



(42) 



to leading order in {Sm2i/ 611112) ^^'^ ^3- ^.s expected, 
the increment R scales with the square of the neutrino 
energy and the VEP breaking scale \(f>\6j which may 
be positive or negative . Since we have restricted \4>Sj\ 
to be much smaller than the atmospheric oscillation scale 
[Eq. (|25|l ]. 101^7 interferes substantially only with the so- 
lar neutrino mass splitting 5m\i, and hence only Sm^i 
appears in the leading order term of R. Naturally, the 
degree of Ve i^t VEP mixing, sin 2(/?, is also a con- 
tributing factor. 

The quantity B depends on the mixing angles of Um 
(some of which we assume can be established from exist- 
ing experimental data), as well as the intrinsic CP phase 
6 and one of the relative phases f3. In the simplest case 
where P = 0, B becomes independent of both 6*3 and i5 
such that 



APt 
A^ 



AE^ sin 6*2 
Sm2i sin 2di 



|0|(57 sin2(/?, 



(43) 



to zeroth order in {Sm2i/5m'^2)- Figure shows the ap- 
proximate equation plotted as a function of \(t>\Sj 
with the neutrino energy fixed at 13 GeV for maximal 
VEP mixing, sin2(/3 = 1. Also displayed in juxtaposi- 
tion are four sets of data points corresponding to exact 
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FIG. 3: AP / AP versus \<j>\5'y for different neutrino energies. 
Maximal Ve i^t VEP mixing is assumed, and (3 — 0. The 
top (bottom) set of curves is for negative (positive) |(;6|(57's. 
We use the best-fit standard oscillation parameters and S3 = 
0.113. 
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FIG. 5: T-trajectory diagram for maximal ^ i't VEP 
mixing, with /3 = 0. The ellipses labelled V + (— ) are the tra- 
jectories in vacuum for a positive (negative) Sm32- Those la- 
belled T+ and T— are the corresponding trajectories in mat- 
ter. Solid (dashed) lines are for a positive (negative) |(^|(57, 
while dotted lines are for the standard, no VEP case. 
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FIG. 4: AP/AP versus |(^|57 for different values of sm2ip, 
while sin'^ ip and (3 are set to zero, and the neutrino energy 
is fixed at 13 GeV. The top (bottom) set of curves is for 
negative (positive) |(/f>|(S7's. 



values of APt/APt calculated numerically for different 
values of S3, 6, and baselines L. Figure |21 is similar to 
Fig. n save for a different neutrino energy, E = 25 GeV. 
Concordance between the exact and the approximate re- 
sults generally obtains. Deviations from the exact results 
appear for the larger |0i57|'s (^ 10~^^) as one tunes up 



the neutrino energy, since satisfaction of the requirement 
H25|l becomes increasingly poor. 

Depending on the relative sign of \(l)\Sj and (Jm^^, the 
ratio APt/APt splits into two branches for any given 
set of E and sin2(^. Since the sign of Sm^i is fixed by 
the solar neutrino data to be positive, the dependence 
is on the sign of |<^|i57 only, as is evident in Figs. ^ and 
12 The top branch, where \4>\Sj < 0, always leads to 
an enhancement in the observed T-asymmetry. The bot- 
tom branch, on the other hand, shows a APt/ APt that 
tends to zero, flips sign, and then continues to grow in 
magnitude as we increase |0(57|. 

More APt /APt versus \(l)\6^ curves are displayed in 
Figs. 13 and ^ for different neutrino energies and mix- 
ing parameters sin 2ip. Figure |31 shows that deviations 
of APt from APt are more prominent for higher en- 
ergy neutrinos, which is consistent with expectation. De- 
pending on the neutrino energy, the enhancement in the 
magnitude of the T-asymmetry in the parameter region 
10^26 < i^j^i < can be anywhere from ~ 20 % to 

a factor of a few. 

It is interesting to note that the corresponding changes 
in the actual oscillation probabilities are generally no 
more than a few percent, as illustrated, for example, by 
the T-trajectory diagrams in Fig. Thus, by compar- 
ing the appearance rates of two T-conjugate channels. 



^ CP - and T- trajectory diagrams were first introduced in Refs. 
I34I and I35I respectively. 
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FIG. 6: /S.Pt / /^Pt versus /3 for different values of 5 for max- 
imal Ve «-> i't VEP mixing. The neutrino energy is fixed at 
13 GeV and |0|(57 = 1 x 10"^^ 
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FIG. 7: Same as Fig.|3 but for maximal Ue ^ VEP mixing 
and Q = 0. 



we will be able to probe a whole new region of VEP pa- 
rameter space that is inaccessible to conventional single 
channel appearance or disappearance experiments run at 
the same energies. 

Before concluding this section, let us remark on the rel- 
ative phase f}. Figure El shows APr/ APt as a function of 
/3, according to Eqs. (|1T|| and for several values of 
(5's at a fixed i<^^|0|(57sin2(p. The message here is clear: 
extra phases in the Hamiltonian do not always contribute 
to increasing the extent of T- violation, particularly when 
6 is already close to the maximally T-violating value of 



: 13 GeV L = 3000 km 




I4)5yI = 1 X 10"- 



sin2tp - 1 sin \|/ - 



0.005 



0.01 



0.015 



P(Ve 



FIG. 8: T-trajectory diagram for maximal ^ Vr VEP 
mixing with fi as the phase variable. The standard CP phase 
5 is set to identically zero. The ellipses labelled V -I- (— ) 
are the trajectories in vacuum for a positive (negative) 5m32 . 
Those labelled T-\- and T— are the corresponding trajectories 
in matter. Solid (dashed) lines are for a positive (negative) 
101(57. For the purpose of comparison, we have also included 
the standard T-trajectories (as functions of S). These are 
represented by the dotted lines. 



7r/2 or 37r/2. Interference between the phases means that 
the final T-asymmetry may even be suppressed to zero. 
On the other hand, if |sin(5| <C |sin/3|, the latter be- 
comes the dominant phase, and the T-asymmetry gener- 
ated therefrom peaks at /? ~ 7r/2,37r/2, as is evident in 
Fig. El We will revisit this issue in a later section. 



B. Pure 



VEP mixing 



This case corresponds to setting V' ~ tt/2. The T- 
asymmetry ratio is given by 



1 



\4>\5^ sin 2(/3, 



where 



A 



C2 sin((5 — a) 



SiCi 



sin (5 



1 



'21 



'32 



52 Sm^i sin a 

53 6171^2 sin (5 ' 



(44) 



(45) 



to leading order in {5m\i / 5111^2) a-nd S3. Comparing Eqs. 
(|^ and we see that since S2 — C2 — 1/V2, the ef- 
fects of this VEP mixing channel on the T-asymmetry 
must be very similar to those from pure Vf. <-> VEP 
mixing. In the case of a = /5 = 0, the changes induced 
in APt are identical for the two mixing channels upon 
replacing \4)\5^ — > — 1(/)|(57 in Eq. 144|) . The actual oscilla- 
tion probabilities, however, are quite different. Compare 
the T-trajectories in Fig. [3 with those in Fig. O 
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Pure ^ Vt VEP mixing 



The case of 
Lp — such that 



mixing is equivalent to setting 



\(f)\Sj sin 2-0, 



(46) 



where 



1 (5to21 *2 sin(Q; — /3 — (5) — sin(a — P + S) 
S2C2 6171^2 sin (5 

S3 sin(Q! — P) 



+ - 



SiCi 



sin (5 



(47) 



is suppressed by a factor of {Sm2i/ 6171^2) or ^3- This im- 
plies that the magnitude of the increment R is generally 
some 0(10^^) — > 0(10^^) times smaller than would be 
observed in the case of, say, pure <-> VEP mix- 
ing for the same \4>\S"/ and mixing parameter, and at 
the same energy. At E = 13 GeV and \(j>S'^\ ^ 10"^^, 
the change in the T-asymmetry is no more than one per 
cent, unless | smS\ <C | sin(a — Thus T- violation ex- 
periments run at energies of a few tens GeV are unable 
to probe VEP in the {v^,Vt) sector to below '-^ 10"^^, 
let alone supersede the \4>5^\ < 10^^^ bound imposed by 
the Super-Kamiokande atmospheric neutrino data |l2j| . 



V. T- VIOLATION WITH 5 = 

An interesting consequence of a nonzero a and/or (3 
is that either or both of these relative phases alone can 
give rise to T- violation in the oscillation probabilities, in 
the absence of any explicit CP- or T- violating phases in 
the mixing matrices U'^j and U'q. Consider for exam- 
ple VEP mixing only in the (z^ei^^r) sector, i.e., "0 = 0. 
From Eqs. (|41|) and (|42|l . the asymmetry APt is given 
approximately by 



APt — Q sin /3x I6S1C1S2C2S3C3 



1 (5m^2'^"^23'^™31 



A12L 



^12Lj , L-i.23 

X sm — - — sm — — 



8E^ A12A23A31 



. A31L 
sm — - — 



(48) 



where 



Q 



2E^\(f>\5j sin2(^ 



^21 



S2 



SlCi 



(5m^2 



617121 \ C2 5m2i"' 



S3 Stti^ 



32 



in the limit S 
T-asymmetry, 



(49) 

0. Comparing this with the standard 



APt — sin (5 x 16 S1C1S2C2S3C3 



1 



'12'' 



'23'^ 



'31 



8E^ A12A23A31 



. A12L . A23P . A31L 

X sm sm sm 



(50) 
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FIG. 9: survival probability at various energies and \<j}5"f\ = 
1 X 10"^^ Solid (dashed) lines are for maximal *-* Vr VEP 
mixing with a positive (negative) |0|(57. Dotted lines are for 
the standard, no VEP case. 



we see that the relative phase /3 can mimic the role of the 
standard phase (5. 

Assuming (5 = 0, Fig. |H1 shows the T-trajectories for 
this case at £; = 13 GeV and \(\)6^\ = 1 x 10"^^, with 
/3 as the phase variable. Because of the factor Q, the 
shape of these T-trajectories can be quite different from 
the standard T-trajectories (as functions of (5), which are 
displayed also in the same figure. 



VI. COMPATIBILITY WITH THE 
SUPER-KAMIOKANDE ATMOSPHERIC 
NEUTRINO DATA 

Since the sensitivity of an experiment to VEP-LIV in- 
creases with T-P, it is instructive to reexamine our case in 
the light of the Super-Kamiokande atmospheric neutrino 
data, whose parent neutrino distribution spans some four 
decades in neutrino energy, P 0.1 — > 10^ GeV. We do 
not have the technology to perform a full-fledged three- 
flavour fit to the data to ensure that the VEP parameters 
considered in this work are not in conflict with experi- 
mental data. However, the same purpose may be more 
than amply achieved by requiring the new survival 
probabilities not to deviate too much from that given by 
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FIG. 10: Same as Fig. El but for |057| = 5 x 10"^'' and 
maximal ^ i/^ VEP mixing. 



FIG. 11: Same as Fig. but for \(t>5'y\ = 2 x 10"^'' and 
maximal ^ Vr VEP mixing. According to Ref. p^ . the 
solid and dashed curves are allowed by the Super-Kamiokande 
atmospheric neutrino data. 



the best fit parameters. 



VII. CONCLUSION 



Figures El and ^1 show the survival probability 
at several energies, with maximal VEP mixing in the 
(v^tI't) and (z/gji'/i) sectors respectively. In the case of 
i/g VEP mixing, there is virtually no deviation from 

the standard survival probability &t E = 10, 100, 
and 500 GeV for \(l)6-f\ = 1 x lO^^s. The case of 
^ i^fi VEP mixing seems somewhat more restricted. 
At E = 500 GeV and \(j}S-/\ = 5 x IQ-^e, the departure 
from the standard probability can be as large as ~ 15 %. 
However, bearing in mind that the neutrino source peaks 
at ^ 100 GeV, this divergence at 500 GeV is most likely 
still acceptable. The case of VEP mixing in the (i^^, i/r) 
sector was considered in Ref. ^l^ in a two-flavour analy- 
sis. For the purpose of comparison, we plot in Fig.llllthe 
t'p survival probability for a just above the upper 

limit of 10^^^ from the said analysis. Clearly, the curves 
in Figs. 1^ and IIUI are no worse than these "acceptable" 
ones. Thus we can be quite confident that the VEP pa- 
rameters considered in this work have not already been 
ruled out by Super-Kamiokande. A proper analysis is 
required to settle this issue. 



We have considered in this work the prospects for ex- 
tending the limits on equivalence principle and Lorcntz 
invariance violations (VEP and LIV) in future neutrino 
factory experiments. We have examined three neutrino 
flavour oscillations, treating VEP-LIV as a sub-leading 
mechanism to the standard mass-mixing mechanism, un- 
der the simplifying assumption that the neutrino-gravity 
couplings of two of the three neutrino gravitational eigen- 
states are nearly the same. The presence of nondegener- 
ate neutrino-gravity couplings introduces additional mix- 
ing and phases which can lead to CP- and T-violation 
even if the phase S in the MNS mixing matrix van- 
ishes. T-violation measurements have the potential to 
yield the most stringent bounds on VEP-LIV. Our case 
study shows that, with suitable neutrino path lengths 
of several thousand kilometres, the limits on VEP-LIV 
for the (I'eji'fj.) and (i^e,i/^) sectors can be lowered to 
\(j)Sj\ ^ 10~^^, which matches the limit for the 
sector imposed by the Super-Kamiokande atmospheric 
neutrino data. At present, we see little prospect for ex- 
tending this limit much beyond the 10~^^ level. 
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